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Heavy-Light Meson Semileptonic Decays with Staggered Light Quarks 
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We report on exploratory studies of heavy-light meson semileptonic decays using AsqTad light quarks, NRQCD 
heavy quarks and Symanzik improved glue on coarse quenched lattices. Oscillatory contributions to three-point 
correlators coming from the staggered light quarks are found to be handled well by Bayesian fitting methods. 
B meson decays to both the Goldstone pion and to one of the point-split non-Goldstone pions are investigated. 
One-loop perturbative matching of NRQCD/ AsqTad heavy- light currents is incorporated. 



1. Introduction 

A major systematic error in many studies of 
heavy-light meson decays has been the uncer- 
tainty due to chiral extrapolations in the light 
valence and light dynamical quark masses. These 
chiral extrapolation errors can be reduced sig- 
nificantly by working with light quark actions 
that have good chiral properties and this has 
motivated us to initiate a program to study 
heavy-light physics using improved staggered (or 
equivalently improved naive) light quarks coupled 
with NRQCD heavy quarks [1]. Our first re- 
sults on B and D meson decay constants on the 
MILC dynamical configurations were presented 
by M.Wingate at this conference [2]. Here we 
discuss initial semileptonic decay results. 

Our semileptonic studies are at a much more 
preliminary stage compared to the decay con- 
stant calculations. We report on quenched sim- 
ulations with light quark mass around Tn s t ran q e 
and the heavy quark mass around the 6-quark 
mass. We employ the Asqtad light quark ac- 
tion, the 0(1/M) and 0(a 2 ) improved NRQCD 
heavy quark action and the tree-level Symanzik 
improved glue action. We work with 200 12 3 x 32 
quenched configurations with a -1 « 1.0 GeV. 



2. Fitting Oscillatory Three-point Corre- 
lators 

Two- and three-point correlators involving 
naive or staggered propagators have contributions 
that oscillate with time and this makes it more 
challenging to carry out fits and extract ground- 
state contributions. Experience gained recently 
with NRQCD/staggered heavy-light two-point 
correlators have demonstrated the usefulness of 
Bayesian fitting methods [1,3]. Hence, much of 
the effort in this exploratory study of semilep- 
tonic decays has focused on testing Bayesian fits 
to oscillatory three-point correlators. The basic 
correlator of interest is, 

C (3) (p.,p s ,t,T B ) = 

££(<M0)J^)$W,T S )) 



* v 



(1) 



In the present simulations the B meson decays 
at rest (ps = 0) and calculations have been per- 
formed with Tb = 12 and Tb = 16. The location 
of the current insertion t was varied between 1 
and Tb- The three-point correlator is fit to the 
form, 
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C^(p v ,p B ,t,T B ) -> 

JV„-1 JVb-1 

^2 (-l) fe * ( * _1) (-1)-J'*( T b-*) 
fe=0 j=0 



0)/ 



V0PS_J0 (Oil) 

NEXP =1+6 
T B = 12 
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In this ansatz every second exponential comes 
with an oscillatory amplitude in both the pion 
and the B channels, We were able to get good 
fits with N v = 1 or 2 and N B = 3 to 8. Fur- 
thermore we found it useful to do simultane- 
ous fits to and to the B meson two-point 
correlator. Fits were done to all data between 
t = 2 and t = T B — 1, i.e. to all but the 
first and the last timeslices, those timeslices on 
which the pion and the B meson respectively were 
created. Figs. I & 2 give examples of fits to ma- 
trix elements of the temporal (Vq) component of 
the heavy-light vector current. The plots show, 
CWfa,t,T B ) e^*- 1 ) e ^ 0) (TB-t). By factor- 
ing out the groundstate contributions, the op- 
posite parity oscillatory terms and other excited 
state contributions show up more clearly. Despite 
the presence of these excitations, the Bayesian fit- 
ting methods allow us to extract the groundstate 
amplitude of interest, A 00 in eq.(2). 

The two figures are for the cases when the 
pion has momentum (0, 1, 1) in units of 2n/(aL). 
Good x 2 /dof were found with N v = 1 and 
N B = 6. The difference between the T B = 12 
and T B = 16 results is typical for all matrix el- 
ements. With T B = 12, statistical errors are re- 
duced and more accurate fits can be made. On 
the other hand, with T B = 16 a hint of a plateau 
can be observed (which will likely disappear with 
higher statistics and a more expanded vertical 
scale). Since with multi-exponential Bayesian fits 
the presence of a plateau is not required, one is 
better off with the more accurate T B = 12 data. 
Nevertheless, working with two different T B val- 
ues provides a good consistency check, which we 
exploit in the next section. 
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Figure 1. Fit results for the groundstate am- 
plitude A o from the (Vb) three-point correla- 
tor with T B = 12. The pion momentum is 
(0, 1, l)27r/(aL). N exp gives N n + N B . 
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Figure 2. Same as Fig.l but for T B = 16. 



3 



3. Results for Form Factors 



Given the amplitudes ^4oo from fits to three- 
point correlators one can determine various form 
factors. In these preliminary investigations we 
find it most convenient to work with the form 
factors /|| and f± defined as, 



/II 
f± 



A OQ (Vo) 
Aoo(Vfe) 



v 



d 01 



\J Zy k 



(3) 
(4) 



and denote the groundstate amplitudes 
of the pion and B correlators respectively. The 
Zy^ are the matching coefficients for the heavy- 
light currents. We use one- loop perturbative re- 
sults for them. 

In Fig. 3 we show results for f\\ and f±. Our 
main results are the octagons and diamonds, 
which correspond to = 12 and Tb = 16 re- 
sults respectively, for decays into local Goldstone 
pions. One sees consistency between the two Tb 
results, with, as expected, smaller statistical (and 
fitting) errors when Tb = 12. 

We have also considered B decays to one of the 
point-split non-Goldstone pions, namely the one- 
link pion. The masses of the one-link and local 
pions differ as M 7r (one — link) / {local) — 1.13. 
The "bursts" in Fig. 3 show form factors for the 
decay to the one-link pion. One sees that at low 
momentum agreement is reasonable with the lo- 
cal pion results. Things start to differ at higher 
pion momentum, however the data has started to 
deteriorate at those momenta and one also ex- 
pects lattice artifacts to be worse there. Things 
should improve considerably as one goes to lat- 
tices finer than the a -1 « lGeV lattices of the 
present study. One can conclude from Fig. 3 that 
taste-breaking lattice artifact errors are smaller 
in form factors than in the corresponding pion 
masses, particularly at low pion momenta. 

The above form factor results come from the 
lowest order (in 1/M) heavy-light currents. We 
have also analysed matrix elements of all the di- 
mension 4 current corrections that come in at the 
next order (two additional currents for Vq and 
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Figure 3. Semileptonic decay form factors /|| and 
f± versus the pion energy E„ in the B rest frame. 



four more for 14). We find that their contribu- 
tions arc small for B mesons giving at most a 
~ 2 — 3% correction. 

4. Summary 

Exploratory quenched studies of heavy-light 
meson semileptonic decays with Asqtad light 
quarks look promising. Simulations are underway 
to calculate form factors on the MILC dynamical 
configurations. 
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